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Abstract 

■  We  present  a  systematic  procedure  for  the  direct  calculation  of  the  free 
energy  and  Its  first  and  second  derivatives  for  the  Islng  model  in  one  to 
three  dimensions  with  a  wide  class  of  synmetry  properties.  This  renormalization 
group  method  is  based  on  the  modified  Kadanoff's  variational  method  (MKUM). 

This  method  is  not  only  general  but  also  very  accurate  numerically  both  near 
and  far  from  the  critical  point.  Further  it  includes  correction  to  scaling 
effects  not  present  in  the  standard  linearized  renormalization  group  treatment. 
This  work  describes  the  techniques  and  presents  some  illustrative  results  for 
the  square  lattice  and  body-centered  cubic  lattice  with  nearest  neighbor  inter¬ 
actions.  A  full  treatment  of  our  results  will  be  given  elsewhere. 
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1 .  Introduction 

In  recent  years,  the  position  space  renormalization  group  (PSRG)  method 
and  its  application  to  phase  transitions  of  various  spin  models  has  been 
widely  studied.  One  important  example  is  Kadanoff's  lower  bound  renormaliza¬ 
tion  group  transformation  (LBRGT)  [1]. 

In  Kadanoff's  LBRGT,  one  has  variational  parameters  in  the  renormaliza¬ 
tion  group  transformation  (RGT)  equation  relating  transformed  and  original 
unit  cell  potentials.  The  optimum  variational  parameters  must  be  chosen  in 
order  to  get  the  best  lower  bound  free  energy  from  the  transformed  Hamiltonian. 
Several  different  methods  with  different  degrees  of  complexity  [2,3,4]  have 
been  proposed  to  determine  these  variational  parameters.  One  such  method  is 
the  Modified  Kadanoff's  Variational  Method  (MKVM)  [3]. 

In  the  MKVM,  the  variational  parameter  p  is  determined  by  minimizing  the 
single  cell  free  energy.  This  leads  to  an  analytic  nonlinear  equation  for  p. 

In  previous  papers  [3,5],  it  was  found  that  the  MKVM  is  very  accurate  for  two¬ 
valued  spin  models  without  external  magnetic  field  in  one  or  two  dimensions.  It 
has  also  been  applied  to  the  three  dimensional  Ising  model  on  a  BCC  lattice 
where  the  derivatives  of  free  energy  were  calculated  by  numerical  differentia¬ 
tion  of  the  free  energy  [6]. 

In  this  paper,  we  formulate  a  systematic  procedure  to  directly  calculate 
the  free  energy  and  its  first  and  second  derivatives  with  respect  to  temperature 
and  external  magn'tic  fields  for  a  wide  class  of  Ising  models  based  on  the  MKVM. 
In  this  way,  we  avoid  the  numerical  errors  arising  from  numerical  differentia¬ 
tion  and  thus  have  more  accurate  derivatives  of  the  free  energy  for  further 
analysis . 
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In  the  application  of  R6  theory  to  phase  transitions,  one  usually  formu¬ 
lates  the  RGT  equations  for  the  coupling  constants,  then  solves  them  for  the 
fixed  point  solution(s)  and  finally  expands  the  RGT  equations  around  the  fixed 
points  to  find  the  linearized  RGT  equation  and  Its  eigenvalues  and  hence  critical 
exponents.  Instead  of  that,  we  apply  the  RGT  to  the  Hamiltonian  of  the  given 
system  and  calculate  the  free  energy  and  its  first  and  second  derivatives  directly. 
From  the  behavior  of  the  variational  parameters  and  coupling  constants  (which 
tend  to  larger  or  smaller  values),  we  can  determine  the  critical  point  of  the 
given  system,  which  has  more  direct  physical  significance  than  the  fixed  points. 

By  analyzing  the  values  of  the  first  and  second  derivatives  of  the  free  energy 
near  the  critical  point,  we  can  determine  critical  exponents,  critical  amplitudes 
and  other  critical  parameters  including  their  less  divergent  corrections  which; are 
of  considerable  current  interest  [7,8].  This  analysis  will  be  presented  elsewhere. 

Since  the  quantities  of  direct  interest  here  are  feee  energies  and  their 
derivatives,  instead  of  the  fixed  point  and  linearized  RGT,  it  is  not  necessary 
to  carry  out  the  calculation  in  coupling  constant  space.  In  fact,  to  write  a 
single  computer  program  applicable  to  many  systems  with  different  symmetry 
properties,  it  is  more  convenient  to  carry  out  the  calculation  in  cell  potential 
space.  This  is  related  to  coupling  constant  space  by  a  linear  transformation. 

The  present  paper  is  organized  as  follows:  In  Section  II,  we  derive 
formulas  for  free  energies  and  their  first  and  second  derivatives  with  respect 
to  temperature  and  external  magnetic  field.  These  formulas  are  written  in  cell 
potential  space  instead  of  coupling  constant  space.  In  Section  III,  we  con¬ 
sider  some  typical  systems  with  various  symmetry  properties  In  1  to  3  space 
dimensions.  We  show  that  we  may  use  very. simple  criteria  to  classify  all 
possible  spin  configurations  on  a  unit  cell  into  different  degeneracy  groups, 
each  of  which  corresponds  to  an  independent  cell  potential.  The  systems  con¬ 
sidered  are  listed  in  Table  I. 


We  illustrate  our  formalism  in  section  IV  by  calculating  the  free 
energy  and  its  first  and  second  derivatives  with  respect  to  temperature  and 
external  magnetic  field  for  Ising  models  on  two  dimensional  square  and 
three  dimensional  BCC  lattices.  Our  results  are  compared  with  exact  [9,10] 
and  series  expansion  values  [11,12]  and  are  shown  in  Fig.  1  and  2.  In 
section  V,  we  give  a  brief  discussion  of  possible  extensions  our  method. 
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11 .  General  Formalism 

Let  us  consider  N  Ising  spins  o^  (=+l)  on  a  general  d  dimensional 
lattice,  which  interact  with  each  other  via  a  Hamiltonian  of  the  form  [1]: 

pH >  (1) 

where  r  is  a  sum  over  d  dimensional  hypercube  unit  cells  of  the  lattice. 

ft' 

°R'  =  ^al' — az^  z s ^  are  sP^ns  corners  of  the  unit  cell  R'  and 

u  (oR>)  is  the  interaction  potential  of  the  spins  crR^  within  the  unit  cell  R'. 

After  a  Kadanoff's  LBRGT,  the  transformed  Hamiltonian  for  N**  Ising  spins 
y.  (=+l)  on  lattice  with  double  lattice  spacing  has  the  same  form  as  Eq.  (1), 
i .  e. : 

(2) 

where  N"  is  the  number  of  spins  on  the  new  lattice  (N'=N/Z),  uR^  (uR)  is  the 
interaction  potential  for  the  Ising  spin  yR=^u^ > — Vz^on  a  hypercube  unit  cell 
R  with  doubled  the  original  lattice  spacing.,  and  uR-  is  related  to  uR(oR^)  by 
the  RGT  equation  (hereafter  we  will  drop  the  subscripts  on  uR-(aR^)  and  uR(yR): 

ocp['v'o4)‘)  =  £  ■e*f[zviO+pzjUi<rI  ~  p)}  /  (3) 

with  u  (o,p)  given  by:  ^ 

uc«r,  p>=  er*‘  +  £ps‘ )  »  (4) 

where  S-jJ:  and  p  is  a  parameter.  It  follows  from  Eq.  (3)  and  (4)  that 
i=l 

u'(u)  and  u(o)  have  the  same  symmetry  properties  with  respect  to  the  point 
group  transformation  of  spins  within  the  hypercube  unit  cell.  Thus  u  and 
u'  can  be  expressed  in  terms  of  the  same  invariant  functions  of  Ising  spins 
on  the  hypercube  unit  cell: 


V(c0=i?-cK;  ’ 


V 


where  g0=(o)=go(y)=l .  The  other  are  defined  and  discussed  in  Appendix  A. 
Let  £={K0,Kj, — Kp  and  £'=(K',Kp — Kp  denote  vectors  of  the  coupling 
constants.  Then  Eq.  (3)  may  be  considered  as  a  transformation  of  the  coupling 
constants:  v 

K'=F  L*,?)  . 

For  Z  Ising  spins  on  a  hypercube  unit  cell,  there  are  2Z  possible  spin 
configurations  and  thus  2Z  possible  hypercube  cell  potential  u'(y).  However, 
many  u'(y)  are  equal  due  to  the  symmetries  of  the  spins  on  a  unit  cell.  In 
fact,  the  number  of  independent  u"{y)  is  the  same  as  the  number  of  independent 
coupling  constants  in  Eq.  (5),  ie  £+1.  Let  denote  the  vector  of  these 
independent  cell  potentials  (u^u^, — ’) .  Then  v'  is  related  to 
K'=(K',Kp — Kp  by  the  linear  transformation  of  the  (£+1)  x  (£+1)  matrix  T: 

W.  )=  ( Xo,  k,'/- .x('  )T  , 

or  briefly  denoted  as: 

V'  =  k'T  , 

T  can  be  calculated  from  Eq.  (5). 

We  also  have  the  inverse  transformation 


where  T-^  T=I. 


k'  =  T/'T 


Now  the  free  energy  calculated  from  H"* (yi-j  is  always  a  lower  bound 

for  free  energy  calculated  from  H(o^  .Ogt—op  and  we  must  vary  p  to  obtain 
the  optimum  lower  bound  free  energy.  In  the  MKVM,  the  variational  parameter 
p  is  determined  by  minimizing  the  single  cell  free  energy.  This  leads  to 
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which  is  solved  at  each  iterative  step  to  determine  P.  From  Eq.  (3),  (4) 

and  (10),  it  is  easy  tt^how  that  Eq.  (10)  may  be  rewritten  as: 

V?r  C 2 Uv  )Tt>  ^ •  P )  _  D  (111 

—  - - c -  -  °  > 

With  F(S^,^,p)  given  by  l/<)  F  ,  f) 

2  F($i,i;  p)=  Cztoskrf/ztosMps^i)  ,  (12) 

Where  S,,,  is  I  oj  evaluated  at  the  i  th  degeneracy  group.  Di  and  ui  are 
1  1  j=l 

the  number  of  elements  and  cell  potential  for  the  i  th  degeneracy  group  respec¬ 
tively. 

To  begin  with,  we  use  (a^  )  of  Eq.  (1)  as  u(a)  in  the  right  hand  side 
of  Eq.  (3)  with  p  given  by  Eq.  (11)  and  calculate  the  cell  potential  u'(u) 
from  Eq.  (3).  This  constitutes  the  first  step  of  the  RG  transformation.  We 
then  use  u"(p)  thus  obtained  as  input  in  the  same  procedure  to  calculate  the 
transformed  potential  u"(p).  This  RG  transformation  is  iterated  further  so 
that  a  series  of  cell  potentials  v'(u),  u"(p), — u ^  (p)  (and  hence  coupling 
constants  K',  K" , — and  the  corresponding  variational  parameters 
p',p", — p^)are  obtained.  In  this  step  by  step  RG  transformation  (a) 
for  v>|  and  p^  tend  to  diminish  to  zero  with  increasing  a  when  the  system 
is  above  the  critical  temperature  and  tend  to  grow  when  the  system  is  below 
the  critical  temperature.  Thus,  after  a  large  number,  say  a,  of  RG  transforma¬ 
tions,  £  times  the  free  energy  per  spin  for  the  original  lattice  may  be 
approximated  by: 


^<f'>  st-  jx  (  k£°+  KZ  ) 
T7Tc  (or/3<(?t,  )  and 


for 
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1 


for  T  <  T£  (or  8  >  Bc)»  where”o0  is  the  spin  configuration  at  T=o. 

To  calculate  the  internal  energy  and  spontaneous  magnetization  of  the 
system,  we  must  take  the  first  derivative  of  f^  with  respect  to  8  and 
h(=8B)  respectively.  Let  q  denote  either  6  or  h.  By  the  chain  rule,  it  is 
easy  to  see  from  Eq.  (13)  that: 

i£!_  li'-il!’,  .  ,  iltL.  .  .  -*3  ?  ,  (14} 

l.lh  A  W 


where 


c  =  ci,  o,-  -  • 


for  T  >  T  and 
c 


c  -  C ! ,  ^  i { 6© )  /  c)2  (  <0 ,  > 


(15a) 


(15b) 


for  T  <  T  . 
c 


it  will  become  clear  later  that  in  order  to  deal  with  systems  with 
different  symmetry  properties  in  a  unified  way,  it  is  convenient  to  carry 
out  the  calculation  in  the  cell  potential  space  u^a^(y)  instead  of  the  coupling 

constant  space  K^a\  From  Eq.  (8)  and  (9),  we  have 

E?u>  ,  t3<*>  t-H  (16a) 


'> 

U 


From  Eq.  (8)  and  (9),  we  have 
'  I  ^(n,,)  1 

s  t> 


(16b) 


Thus  Eq.  (14)  may  be  rewritten  as: 

lvl *> 

^  t  &  B  o  Vc 


f  (17) 


i 


In  Appendix  B,  we  derive  very  simple  general  formulas  for  T  C  that  allow  us  to 
avoid  explicitly  calculating  the  matrix  T. 

It  should  be  noted  that  for  q=h  in  Eq.  (14),  for  m=l , — ,a  must 
include  all  odd  spin  coupling  constants  even  if  we  evaluate  the  row  vector 

gjp  and  matrix  at  points  where  the  odd  spin  coupling  constants  vanish 


These  points  are  equivalent  to  the  points  in  the  cell  potential  space  with 
u  (o- ,.-..oz)  =  u  (-c^,..-o2).  Thus  to  evaluate  the  spontaneous  magnetization 
in  cell  potential  space,  we  must  consider  u(o^,  °2 anc*  u  (-ai»  ~°2  ”~a 
as  independent  cell  potentials. 

To  calculate  the  specific  heat  and  susceptibility  of  the  system,  we  must 
take  the  second  derivative  of  f^  with  respect  to  6  and  h.  From  Eq.  (17), 


it  follows  that 


py*’.  r  V’ %&>  .  ,  tv1""  . 

st1  1  U3-  tv*  TF7” 


+  i£" 

“t' 

-t  -  -  * 


^v1 


4-^  *>  VU  > 

b%  ' 


* 

£  S 


t#-1’ 

^col) 


±^.kkL{>. . .  . . .  L 

s:  nt<»>> 


By  the  chain  rule  A.  ■%.  *  for  m=|  ♦••e*  is  further 
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|£J’  and|l 

{>%  Tv- 


1  rU » 


Thus  to  calculate 


,  we  must 


'i 


evaluate 


&tj°  r-v'° 


Si/< 


5 


’  su0^ 


and 


(c>  > 


for  i,  j,  k  =  1,  — ,  4+1.  It  is  straightforward  to  calculate  the  first  two 


quantities  from  the  initial  cell  potential  of  Eq.  (1).  From  Eq.  (3),  we  have 

i-4  u(<r,p>)]  , 

where  (y-j.--uz)  is  one  of  the  spin  configuration  belonging  to  the  j  th 


degeneracy  group  and  u(a,p)  is  given  by  Eq.  (4).  The  variational  parameter 

of  Ec 
(m-1) 


p  of  Eq.  (20)  is  essentially  determined  by  the  initial  cell  potentials  , 

S'U? 

u'  '  '.  This  fact  must  be  taken  into  account  in  the  calculation  of  y^-T). 

0  y  * 

Thus  we  have 

jyr’.  =  ( 1  )  , 

tv;.""’  1  f-?1--"';)  1  3f  jh';'-1’ 

Where 


(Jls£L)  = 


Z  -f  €<p(z  P^VrH^p)] 

Z  eAp[2  pz  l'k4-U(<f,p)] 


1^=  £  **pl>  !/*-&)+  pi  ^-u(^r)]( fr 


2-^ptZ  17U'l)C^)+pZ>k^-M  (cr,  p )7 

<r  ksj 


and 


BP  _  _  Z  +^p(z  Tp 


i+» 


(k-i  1  a* 


(20 


(21 


(22 


(23 


* 


where  F(S-|,..;p)  is  given  by  Eq.  (12).  in  Eq.  (22),  is  a  sum  over  all  con¬ 
figurations  of  a  corresponding  to  the  cell  potential  u^m’^(o).  Eq.  (24) 
is  derived  from  Eq.  (11)  and  the  notation  is  the  same  as  that  of  Eq.  (11). 


From  Eq.  (21),  it  follows  that 


-(i 

£  ofj  _  /  h  _ _ 

^  2)  3  ;p  T  a'li^ap 

,  12^  ^  P  ,  f  W?' 

'+  ap 

4  9  p—  ' ]-2£~  . 

.  a  P2-  ’dV^  ^p  ap  aD[-M)J  a,ii(|cM) 

It  is  straightforward  to  derive  equations  for  the  partial  derivatives  on  the 
right  hand  side  of  Eq.  (25).  They  are  very  involved  so  we  do  not  report  them 
here.  It  is  obvious  that  Eq.  (25)  is  symmetric  with  respect  to  the  indices 
K  and  i.  We  have  used  this  fact  to  check  our  expressions  for  the  right  hand 
side  of  Eq.  (25). 


(.HI  )r>  t 


tfr  2SL. 


This  completes  our  formulation  of  the  equations  for  the  direct  calculation 


of  the  first  and  second  derivatives  of  the  free  energy. 


In  order  to  carry  out  the  configuration  sum  I'  in  Eq.  (22)  and  also  the 

equations  for  the  partial  derivatives  with  respect  to  v^~^»  on  t^ie 

right  side  of  Eq.  (25),  we  must  establish  the  correspondence  between  spin 

configurations  (o^,  —  az)  and  the  independent  cell  potentials  i>-(oj,  —oz). 

That  is,  we  must  classify  the  2Z  possible  spin  configurations  into  different 
groups,  such  that  all  configurations  of  the  same  group  have  the  same  cell 
potential  u. (a^ , — az),  where  i  runs  from  1  to  J£+l.  This  correspondence 
depends  on  the  space  dimension  and  symmetry  properties  of  the  system.  However, 
we  can  establish  very  simple  general  criteria  for  the  purpose  of  such  con¬ 
figuration  classification. 

In  Appendix  A,  we  list  invariant  functions  for  certain  systems  in  one 
to  three  dimensions.  We  also  list  the  relations  between  invariants  of  a 
given  system.  From  such  relations,  it  is  easy  to  see  that  all  invariants  for 
a  given  system  may  be  expressed  in  terms  of  certain  basic  invariants.  These 
are  also  given  Table  I.  Thus  by  Eq.  (5),  any  cell  potential  u(o-j, — o  ) 
can  be  expressed  as  a  function  of  these  basic  invariants.  So  we  may  use 
these  basic  invariants  as  criteria  to  classify  the  spin  configurations. 

The  generation  of  all  possible  configurations  and  their  classification  into 
different  groups  based  on  their  basic  invariants  may  be  carried  out  simply 
by  computer.  This  scheme  is  briefly  described  in  Appendix  B. 

In  Table  I,  we  also  list  the  values  of  £+1  and  possible  applications  of  the 
considered  models. 


rw 
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Table  1.  Basic  Invariants  and  possible  applications  for  the  considered  models. 
Note  that  the  definition  of  (see  Appendix  A)  may  differ  from  case  to  case. 


Model 

Dimension 

Basic  Invariants 

A+l 

application3 

SU- 

1 

91 

3 

S2 

2 

92 

5 

Isotropic  SQ  lattice  with  nn 
interaction  only. 

S3 

2 

9r 

92 

6 

Isotropic  SQ  lattice  with  nn, 
nnn  and  4  spin  interactions. 

S4 

2 

9r 

g2’  g3 

7 

Anisotropic  SQ  lattice  with  nn, 
nnn,  and  4  spin  interactions. 

S5 

2 

V 

g2 

9 

Triangular  lattice  with  nn 
and  nnn  interactions. 

S6 

3 

gi 

9 

Isotropic  BCC  lattice  with 
nn  interaction  only. 

S7 

3 

9r 

92*  g3 

22 

.  Isotropic  SC  lattice  with  all 
possible  even  spin  interactions 
within  the  primitive  unit  cell. 

S8 

3 

9] » 

V  93»  g4>  g5 

34 

Anisotropic  SC  (tetragonal) 
lattice.  Crossover  from  d=2 
to  d=3. 

S9 

3 

9r 

g2>  93»  g4»  g5 

46 

Isotropic  FCC  lattice  will  all 
possible  even  spin  interactions 
within  the  primitive  unit  cell. 

a.  nn  =  nearest  neighbor,  nnn  =  next  nearest  neighbor. 
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IV.  An  Example:  Permulation  Symmetry 

In  this  section,  we  apply  the  formalism  developed  in  the  previous  sections 
to  an  Ising  model  with  isotropic  nearest  neighbor  interactions  on  1  dimensional, 
2  dimensional  square,  and  3  dimensional  body-centered  cubic  lattices.  The 
Hamiltonian  in  these  cases  may  be  written  as: 


/3  H  =  -  *  £  Si  <r -  h  r<r,-  =  -  2.  , 

HJr3  K<5c  **  dj?)  +  ^  (^ft-  *  (fj  )  +  K<*C  , 


(26a) 


(26b) 


where  R  is  a  unit  cell  with  one  spin  ac  in  the  center,  Z  spins  , — az  at 
the  corners  and  the  factor  0(=  ^)  has  already  been  included  in  K  and  h. 
Performing  a  decimation  calculation  which  sums  over  the  central  spin  in  each 
cell,  we  obtain  an  effective  Hamiltonian  for  the  remaining  spins.  The  result¬ 
ing  cell  potential  is 


Z  (*•■*’“  +  <J£  >  (27) 

Eq.  (27)  has  permutation  symmetry  with  respect  to  a-j ,  o^,  ~~°z  and  Thus  accord¬ 
ing  to  Appendix  A  can  be  classified  as  SI,  S2,  or  S6  for  d  equals  1,  2  or  3 
respectively.  We  can  use  =ct^  +—  +oz  as  the  basic  invariant  to  classify 
the  2Z  possible  configurations  of  ,  --az  into  different  degeneracy  groups, 
each  of  which  corresponds  to  the  same  cell  potential.  We  then  use  R  of 
Eq.  (7)  as  uQ  in  Sec.  2  to  calculate  the  free  energy  per  spin  and  its  first 
and  second  derivatives  with  respect  to  k  and  h.  We  must  divide  all  these 
quantities  by  a  factor  2  because  Eq.  (27)  was  obtained  after  a  decimation 
calculation.  The  results  for  d=l  agree  extremely  well  with  the  exact  results 


obtained  by  transfer  matrix  method  [9].  These  results  also  provide  a  check 
on  our  computer  program. 

A  few  results  for  d=2  and  3  are  shown  in  Fig.  1  and  Fig.  2  respectively. 

Fig.  lb  and  Fig.  lc  should  be  compared  with  Fig.  6  and  7  in  Kadanoff  et  al.'s 
paper  [1].  Their  results  were  obtained  by  numerical  differentiation  for  the  free 
energy  and  thus  involve  more  numerical  error  than  the  present  work. 

In  future  work,  we  will  present  the  results  of  this  method  applied  to 
several  models  with  different  symmetry  properties.  Preliminary  analysis 
indicates  that  we  obtain  very  accurate  values  both  close  to  and  far  from  the 
transition  temperature,  as  is  clear  from  Fig.  1  and  2.  The  application  of 
this  method  near  the  transition  temperature  is  of  particular  interest  since 
it  allows  an  evaluation  of  correction  to  scaling  effects  due  to  nonlinear 
terms  and  all  the  (linearized  renormal ization  group)  eigenvalues. 


V.  Discussion 

From  the  previous  sections,  it  is  clear  that  to  calculate  the  free  energy 
and  its  derivatives  for  given  model.  It  is  sufficient  to  use  the  basic  invariants 
of  such  model  and  not  necessary  to  use  the  relations  among  all  invariants,  i.e. 
the  even  numbered  equations  of  Appendix  A.  For  a  given' system,  to  obtain  the 
former  is  much  easier  than  the  latter.  However,  if  we  desire  to  calculate  the 
flow  of  the  coupling  constants  (k-|,  kg,  "“kg)  In  parameter  space  as  the  step  by 
step  RG  transformation  is  carried  out,  we  may  use  the  even  numbered  equations 
of  Appendix  A  to  calculate  the  U+l)  x  U+l)  matrix  T  and  hence  T_1  of  Eq.  (8) 
and  (9)  respectively. 

The  method  used  in  this  paper  may  be  easily  extended  and  applied  to  systems 
with  symmetry  properties  different  from  those  listed  in  Table  I.  In  particular, 
we  may  combine  the  ideas  of  this  paper  and  our  previous  paper  [14]  to  calculate 
derivatives  of  free  energy  for  anti  ferromagnetic  systems. 
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Figure  Captions 


Fig.  1.  Derivatives  of  the  free  energy  for  the  Ising  Model  on  the  SQ  lattice. 

2 

(a)-f£  as  function  of  K.  (b)  -|^2  as  function  of  K.  (c)  as 
function  of  K,  where  M  =  -|p  (d)  lnX-a$  a  Function  of  k,  where 

_K^h. 

The  solid  curves  in  Fig.  la,  lb  and  1c  represent  Onsager  [9]  or  Yang's 
[10]  exact  solution.  The  solid  curve  in  Fig.  Id  is  obtained  from  high 
temperature  [12]  and  low  temperature  Tn]  serious  expansion  for  K<0.4407 
and  K>0.4407  respectively. 

Fig.  2.  Derivatives  of  the  free  energy  for  the  Ising  model  on  the  BCC  lattice. 

2 

(a)  -|£  as  function  of  K.  (b)  as  function  of  K.  (c)  M  as 

function  of  k,  where  M  =  -|£.  (d)  In  X  as  a  function  of  k,  where 

^2f  * 

“kf2v  The  sollcl  curves  are  obtained  from  high  temperature  [12] 
or  low  temperature  [11]  series  expansion  for  K<0.1574  respectively. 


Fig.  Al.  Location  of  Ising  spin  on  one  dimensional  unit  cell. 


Fig.  A2.  Location  of  Ising  spins  a-j>  on  the  two  dimensional  unit  cell. 

a  and  S  are  the  primitive  translation  vectors  and  a  is  the  angle  between 
them. 


Fig.  A3.  Location  of  Ising  spins  — og  on  the  three  dimensional  unit  cell 

-4. 

a,  6,  and  c  are  the  primitive  translation  vector  and  a,  g,  y  are  the 
angles  between  them. 


4 
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APPENDIX  A 

Symmetry  Properties  and  Invariant  Functions  for  Ising  Spins  on  Hypercube  Unit  Cell. 

The  location  of  the  z  (=2^)  Ising  spins  on  a  hypercube  unit  cell  are  shown 
in  Fig.  Al,  Fig.  A2  and  Fig.  A3  for  d=l,  2  and  3  respectively.  We  shall  use 

K-  to  denote  the  two-spin  coupling  constant  between  o.  and  c.  which  appears 

-  j  1  j 

on  the  right  hand  side  of  Eq.  (5). 

Now  we  consider  the  following  possible  systems  and  write  down  the  invariant 
functions  on  the  right  hand  side  of  Eq.  (5)  for  them. 

51  d=l 

91  =  °l+a2’ 

92  =  al°2‘  (A1) 

They  satisfy  the  relation: 

g22  =  9l2/2  -1  .  (A2) 

52  d=2  u(o^.--a^)  has  permutation  symmetry. 

91  =  °l+a2+a3+a4’ 

g2  =  a2+02a3+a304+0401+a1 °3+a2°4 ’ 

g4  =  a1a2a3o4.  (A3) 

They  satisfy  the  relations 

g2  -  g,2  n  -2, 

9-j  “  ^192  ~  ^9i  3  /3» 

94  =  Cg -| 9 3  -2g2]  /4.  (A4) 

53  d=2  u(o-|,--o4)  has  rotation,  reflection  and  inversion 

symmetries.  K12=K23=K34=K41’  K13=K24 
but  K12  f  K^3  For  the  unit  cell  of  Fig.  A2,(  a{  =  [S|  , 


1  »- 


91  =  °l+02+03+04’ 

g2  =  ffla3+a2a4* 

s3  =  alo2+02a3+a3°4+O4al  ’ 

94  =cfic2a3+a2a3£74+c304al+o4olo2’ 

95  =  ala203°4' 

They  satisfy  the  relations 

93  =  gT%  ‘  92  ~2' 

94  =  9192  -  91 

95  "  92  ^  - 

S4  d=2  u(0-j,--ct4)  has  reflection  and  inversion  symmetries. 

K12=K34’  K23=K41*  K12^23' 

K13=K24*  For  the  unit  cel^  F^9‘  A2,  ®  ^  ^  *  a=90° 

91  =  °l+02+03+04’ 
g2  =  ol03+o2a4, 

g3  =  ol02+o3o4, 

94  3  cl°4  A  °2°3, 

95  =  °l°2a3  +  a2°3°4  +  °3°4°1  +  d4al°2. 

g6  =  0^203^ 

They  satisfy  the  relations: 

942=  gi2/2-g2-g3-2, 

95=9192_1 

'  g6=g22-l 


(A5) 


(A6) 


(A7) 


(A8) 
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S5  d=2  u(a^,--o4)  has  Inversion  symmetry. 

K,2=K23=K34=K41,  K13  t  K24.  for  the  unit  cell  of  Fig.  A2, 
| a |  =  |S|,  a  =  60°. 

92=02+o4* 

g3=0la3, 

94=0204* 

gg=a-j  a2+a2a3+a  304+0401 

96=0la2a3+0l°3a4, 

g7=ala204+02a3a4, 

98=ala2a3°4. 

They  satisfy  the  relations: 

93=9l2/2-l, 

94=922/2-l , 

95=9192’ 

96=9293’ 

97=9]94> 

98=9394 


(A9) 


(A10) 
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56  d=3  u(o^,--Og)  has  permutation  symmetry. 

Let  g.  denote  the  sum  of  all  possible  products  of  i  different 
o's.  for  i=l,.*8  are  the  basic  invariants  of  the  permutation 
It  is  easy  to  show  that  g^’ satisfy  for  relation: 

g^=c^+02+--+0g  , 

Vr^s,  -  (8-1+,)  9i-i]  / ( i+i ) . 

for  i=l,...7.  Thus  each  g..  may  be  expressed  in  terms  of  g-j. 

57  d=3  u(o^,— Og)  of  Eq.  (5)  has  symmetry  properties  of  the 

simple  cubic  unit  cell,  i.e.  0^  point  group. 

g-|=a.j+o2+03+--+0g. 

92=0la2+ —  02  terms), 

g3=a-jO^+ —  (4  terms), 

94~°1°3+ —  (12  terms), 

95~°l°2a3a4+ —  (®  terms), 

gc=010oO.aoQ+---  (24  terms), 


group. 


(All) 

(A12) 


g6=0l°203a8+— 

97=0la203°6+ 


(8  terms), 
(24  terms). 


98=0102°307+ —  (24  terms), 

g9=ala2a7a8+ —  (6  terms)> 

9 10”ala3°6°8+  (2  terms)’ 

9 1  l=°i 02a3°4a5a6+ — ( 1 2  terms)> 
9 12=ala2a3°4a5a7+-~-(l 2  terms)» 
g 13=°la2a3a507°8+ — (4  terms)> 

914r°lO2o304a5a6a7G8(l  term), 
g.  =0,0,0,+—  (24  terms), 


915=°l°2a3+— 
9 16"al°2°7+ 

9  1 7=C1°3°6+ 


(24  terms), 
(8  terms). 


918=°lc^2c3a4c5+—  (24  terms). 


919s°la2°3cW’""”  (24  terms), 

920=ala3°4cr6a8+ —  (8  terms^* 

921=0^0202^02020^+-—  (8  terms). 

Where,  in  order  to  save  space,  only  one  typical  term  for  each  invariant  is 
written.  It  is  easy  to  show  that  the  of  Eq.  (A13)  are  related  to  each  other 
by  the  equations: 

V  ¥/ ' VS3, 

/  ^ 

g6=  g3g4-2g2. 

98=9293-294, 

9?=932/2-2, 

95=922/2"6-294-98‘99, 

97=(9294'492-693-296)/3‘ 

g10~g4^/^~2_394g5^3_98^3”99^3’ 

g14=g102/2"1’ 

grrg2gi4’ 

912=93914’ 

913=94914’ 

g16=9l93'9l, 

915=(gl92"3gr916^/2, 

g17=(g1g4-3grg15-g16)/3. 

918=91491 5* 

919=91491  , 
g20=g14g17, 

921=914gl 


S8  d=3  u{a1 ,  a^,  --Cg)  of  Eq  (5)  has  symmetry  properties  of  the  tetragonal 
unit  cell.  For  the  unit  cell  of  Fig.  A3,  |aj  =  |6|  f  |c( ,  a  =  8  =  y  =  90 

9l  sal  +a2+a3+VWa7+a8, 

g2=al °2+a203+o3a4+04al +05a6+0607+0708+08°l , 

g3=ala5+0206+a307+0408 
g4=°l °7+0208+o3a5+0406  * 

g5=0l 03+o2°4+0507+0608, 

g6=al06+0205+0207+0306+0308+040/0l08+04a5, 

g7=0l °203O4+O506a708 

g8=0l a20605+020307a6+0304a8a7+01 a4a5a8 , 

g9=al020308+ —  (8  terms) 

91 0=al a2°6cr8+ —  (16  terms) 

9n=0l °2°306+ — -  (8  terms) 

g12=al °2a307+ —  (16  terms) 

g13=0la2a6a7+ —  (8  terms) 

g14=al-a207°8+ —  (4  terms) 

g15=0l050307+02a6a408, 

916=0l030608+0204a507, 


917=ol°20304a5a6+— 

(8  terms). 

g18=al°2a3a506a7+— 

(4  terms). 

g19=olc2c3a507a8+-- 

(4  terms). 

g2O=olo2°304a507+— 

(4  terms). 

g2rola20405°6a7+— 

(8  terms), 

g22=0l a20304050607a8 , 

923"0l°2a3+— 

(8  terms), 

g24=0l02°6+— 

(16  terms) 

g25=0r2°7+— 

(16  terms) 

g26=cl  C5C3+ — 

(8  terms). 
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927=01°306+—  (8  terms), 

928=a]°2°3;,,a5+—  (8  terms), 

9  29=01 °2°305a6+  —  (16  terms), 

g3O=0la3040506+ —  (16  terms), 

931=0102030507+ —  (8  terms), 

932=0la203a608+—  (8  terms), 

g33=0l02°304050607+—  (8  terms). 

Where  in  order  to  save  space,  only  one  typical  term  for  some  invariants  is 
written.  Invariants  of  Eq.  (A15)  satisfy  the  relations: 
g6=  g12/2-g2-g3-g4-g5 

9<r  94V2g3 
910=94g6'2g2 

9irg3g5‘2g4 

972=g293“2g6 

913=g394"2g5 

g16=  (g5Z+  g62-4-g13-2g5'942-g2Z-932)/6 
97=(g22+4g}6-g62)/4 
915=95Z/2'2"97"916 
g8=g3Z/2-2"915 
girg//2~z~gi5 
g22~g7  ^2 ~ ^ 

gi =g22qi-  15  for  i=17’-~~21 * 

923=g7  ^92”93-g4+2g5-g6_2^^ 
g24=2g23‘gl(g2+g6'2g5+2g4)/2’ 

925=2923"gl^g2"g6+2g5"2g4^2’ 

926=“2g23fgl  ^g2'g6+2g5"2^2, 
g29=g23"gl  ^g2”g6^2. 


(A1S) 
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9i*9229i-5  for  1  =  28,  29,  •~“32» 

933=92291 

S9  d=3  u(a-],--cJg)  has  the  symmetry  properties  of  the  primitive  unit  cell  for 
the  FCC  lattice.  The  longest  diagonal  is  in  the  ct-j,  a7  direction.  For  the 
unit  cell  of  Fig.  A3,  |a |  =  |S|  =  |c j ,  a  =  6  =  y  =  60°. 

9r°i+o7 

92=02+04+05+o3+o6+08 , 

93=a2a4+a4a5+a2°5+a3a6+a3a8+a6a8, 
g4=0l(03+06+08)  +07  (°2+04+a5), 

95=0208+o3ct5+0406 ,  ' 

96=91°7, 

g7=°l  ( °2+a4+a5 ) +°7 ( a3+°6+a8 ) » 

98=0203+02o6+o5°6+0508+0408+0304 , 


99=0l°203a4+— 

(6  terms). 

9lO=0l°2a3a8+— 

(12  terms). 

9ir0la2°4a7+™ 

(6  terms). 

912=o2a306a8+— 

(6  terms). 

913=0l02°3a6+— 

(6  terms), 

91 4=al 02a4c5+03°6 

Q 

Q 

03 

* 

915=0la203c7+— 

(6  terms) , 

916=°2a3a4a8+— 

(6  terms). 

917=0203°708+— 

(12  terms). 

9 1 8=01 °2a7a8+0i °4 

a6°7+al0503°7 

919_03a4O5c6+a2o3a5o84o2°6a4a8, 
92O=0lo306°8+o2o4o5o7, 
921=0la2°3°4a5a7+ _  ^  terms), 

922=01 C2°3C4C608+ —  (6  terms). 
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g23~0l020305a708  — 

924=o2°30405o6a8, 

(3  terms) 

g25=0l0203°40506+— 

(6  terms) 

g26=0l02o304°607+— 

g27=°l 0203a4a5a607a8 

(6  terms) 

928=al02°4+— 

(6  terms). 

g29=a203°6+— - 

(6  terms). 

g3O=°l°206+— 

(12  terms) 

g31=0l07(02+04+05+03+o 

6+08), 

g32=0i02°8+— 

(6  terms). 

g33=020608+— 

(12  terms) 

g34=0l°306+— 

g35=03a608+02CT405 ’ 

(6  terms). 

g36=0l02a30405+— 

(6  terms) 

g37=0l020304°7+— 

(6  terms) 

g38=ala20304a6+— 

(12  terms) 

g39=02o304°5V— 

(6  terms). 

g4O=0la3a405°6+— 

(6  terms). 

g4r0la304a5a7+— 

(12  terms) 

g42=0i O2°306°8+ — 

(6  terms). 

943=0la20405°7+0l0306a807, 
g44=al020304050608+02030405°60708, 
g45=0i0203°4050607+—  (6  terms). 

Where  in  order  to  save  space,  only  one  typical  term  for  some  invariants  is 
written.  Invariants  of  Eq ( A1 5 )  satisfy  the  relations: 

96*g,2/?-i 

g7=glg2'g4 

g8=g22/2-3'g3'95 
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APPENDIX  B 

Details  of  the  Calculation 

In  this  Appendix,  we  outline  some  techniques  used  in  our  computer 
program  to  obtain  the  numerical  results.  Our  computer  program  can  be  applied 
to  all  systems  listed  in  Appendix  A  with  space  dimension  and  symmetry  properties 
as  input  parameters.  The  important  programming  techniques  are  as  follows: 

1.  Generation  of  all  possible  spin  configurations.  Let  z  denote  the  number 
of  spins  on  a  hypercube  unit  cell  of  the  d  dimensional  lattice  and  NCF  denote 
the  number  of  all  possible  configuration  of  such  z  spins  o^»— az,  then 
Z=2^,  NCF=2Z.  Let  I  be  an  integer  running  from  o  to  NCF-1.  We  express  I 

as  a  z  digit  binary  number  and  let  1  correspond  to  spin  up  (+1)  and  0 
correspond  spin  down  (-1),  then  each  integer  from  0  to  NCF-1  correspond 
to  one  spin  cofiguration  where  integer  0  corresponds  to  all  spin  down 
configuration  and  NCF-1  correspond  to  all  spin  up  configuration  which  is 
the  ground  state  of  the  ferromagnetic  system.  Thus,  with  d  as  input  parameter, 
we  can  generate  all  possible  spin  configurations. 

2.  Classification  of  spin  configurations.  From  Table  A1 ,  it  is  clear  that 
for  space  dimension  d>2  there  are  some  systems  with  different  symmetry 
properties.  Thus,  for  each  spin  configuration  generated  in  the  manner 
described  above,  we  calculate  the  values  of  the  basic  invariants  based 

on  the  parameter  "sym"  and  then  classify  these  spin  configurations  into 
different  degeneracy  groups  such  that  the  configurations  in  the  same 
group  have  the  same  values  for  all  basic  invariants  (and  hence  cell 
potential)  and  in  different  groups  have  at  least  one  different  value.  We 
label  the  degeneracy  group  in  such  a  way  that  the  ground  state  configuration, 
which  is  all  spin  up  configuration  for  ferromagnetic  system,  belongs  to  the 
final  group.  For  each  group,  we  also  store  the  degeneracy  Di  (ie  the  number 
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of  configurations  in  the  group)  and  the  configurations  belonging  to 
this  group.  We  will  use  to  denote  the  vector  whose  elements  are 
degeneracy,  ie  5=(D^,D2,-— D^). 

3.  Calculation  of  the  column  vector  T~^. 

In  this  section,  we  will  derive  formulas  for  the  vector  T  which 
appears  at  the  end  of  Eq.  (17)  and  (18).  From  Eq.  (5b),  it  is  easy 
to  show  that 

K  NCF  i^l  ^  Di* 
o 

Thus,  Eq  (13)  may  be  rewritten  as: 

f(“)=  _  l/z“(l/NCF  u.(aWln2)  (B2a) 

for  T>T  and 
c 

f(“)=  _  i/Za  (  (“j  +  In  2) 

for  T<TC,  where  u^s,+l  is  the  ground  state  cell  potential.  We  now 
take  the  first  derivative  of  Eq  (B2)  with  respect  to  q.  By  the  chain 
rule,  it  is  easy  to  show  that: 


6f^ 

SuO  Sul  .... 

So 

(m) 

1/Z“  t  . 

(B3) 

where 

6q 

6q  *  SuO 

So 

15=1)  " 

U 

v<v 

V' 

■~Vi  >5® 

(B4a 

for 

T>T  and 

c  V<°  • 

0  , 

, - ,  0,  1) 

(B4b 

for 

T<T  .  Comparing 

Eq 

(17)  and  (B3),  it 

is 

obvious 

that  T_1t 

of 

Eq  (17)  is  jsut  tQ  of  Eq  (B3)  and  (B4). 


